epl draft 



Conductance through quantum wires with Levy-type disorder: 
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Abstract. - In this letter we study the conductance G through one-dimensional quantum wires 
with disorder configurations characterized by long-tailed distributions (Levy- type disorder). We 
calculate analytically the conductance distribution which reveals a universal statistics: the dis- 
tribution of conductances is fully determined by the exponent a of the power-law decay of the 
disorder distribution and the average (InG), i.e., all other details of the disorder configurations 
are irrelevant. For < a < 1 we found that the fluctuations of InG are not self- averaging and 
(InG) scales with the length of the system as L a , in contrast to the predictions of the standard 
scaling-theory of localization where In G is a self-averaging quantity and (In G) scales linearly with 
L. Our theoretical results are verified by comparing with numerical simulations of one-dimensional 
disordered wires. 
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j Quantum coherent electronic transport through disor- 
dered conductors has been widely studied from a funda- 
mental and practical point of view. For instance, the pi- 
oneering ideas by Anderson on the localization of elec- 
tion wave functions in disordered conductors and the one- 
parameter scaling approach to this phenomenon [T] have 
been of large impact on condensed matter physics. The 
progress in the theoretical description of the electronic 
transport has been accompanied and stimulated by ad- 
vances in the fabrication of small electronic devices, where 
coherent electronic transport is possible. 

, On the other hand, random processes characterized by 
probability densities (Levy-type distributions) p{x) with 
long tails, i.e., for large x, p(x) ~ l/x 1+a with < a < 2, 
have been found in very different phenomena in nature 
and human activities. For instance, it has been seen that 
the movement patterns of some marine predators follow 
a Levy-type distribution [2]. In economy, fluctuations of 
the stock market indices can be described by Levy mod- 
els [3]. Mathematicians have been investigating the prop- 
erties of this class of probability densities since seminal 
works by Levy [H[5]. Among these properties we remark 
the generalized central limit theorem (GCLT). We recall 
that the central limit theorem states that the normalized 
sum of independent variables with finite mean and vari- 



ance is normally distributed in the limit of a large num- 
ber of variables; the GCLT gets rid of the finite variance 
condition and the resulting limit probability distribution 
is a Levy or a-stable distribution. In physics, Levy-type 
distributions have been applied to several areas such as 
statistical mechanics, fluids, and dynamics of chaotic sys- 
tems [5j . Therefore Levy processes have become of a broad 
interdisciplinary interest. 

Levy processes have been found in different phenom- 
ena, but an experimental study of these processes in a 
controllable way has not been possible until very recently. 
New techniques in the fabrication of materials have al- 
lowed the experimental realization and study of Levy 
transport [51(7]. For instance, in Rcf. [f)j SiC nanowires 
have been fabricated via self-organized processes. It turns 
out that the diameter of the SiC nanowires shows fluctu- 
ations whose distribution follows a power-law decay. In 
Rcf. |7j Barthelemy and collaborators were able to make 
a disordered medium where light travels across glass mi- 
crospheres whose diameter follows a Levy-type distribu- 
tion. Motivated by these kind of experiments, here we 
shall study the quantum coherent transport of electrons 
across one-dimensional (ID) disordered wires with Levy- 
type disorder. Previous theoretical works [5] have pointed 
out and studied the effects of Levy processes in the elec- 



p-1 



F. Falccto V. A. Gopar 



Fig. 1: Schematic of a quantum wire with scatterers (shaded 
regions) randomly placed accordingly to a Levy-type distribu- 
tion. The wire is fed with electrons by the reservoirs (ii t 2 



tronic transport. 

In this work we obtain the distribution of the conduc- 
tance of electrons traveling through a ID disordered wire 
whose configuration of impurities follows a density distri- 
bution with a power-law tail, fig. 1. We assume that the 
electrons travel coherently across the sample. Thus inter- 
ference effects arc reflected in transport quantities, such 
as the conductance, and due to the Levy-type distribu- 
tion of the scatterers novel effects are seen. For instance, 
the coexistence of ballistic and insulating regimes. More 
interesting from a fundamental point of view is that a new 
universality class emerges, as the distribution of the con- 
ductance is completely determined by the power-law tail 
exponent a of the disorder distribution and the ensem- 
ble average (lnG). No further information is needed to 
determined the conductance statistics. 

To study the conductance we adopt the Landauer- 
Buttiker approach in which the electronic transport is seen 
as a scattering problem. In this approach the dimension- 
less conductance G (in units of the conductance quantum 
2e 2 /ft) is given by the transmission coefficient. With re- 
spect to the statistical study of the conductance, we shall 
use some results from the standard scaling theory devel- 
oped in Rcfs. [SHE!]. These works considered the weak 
scattering limit {Xp <C I, being Xp and I the Fermi wave- 
length and the mean free-path, respectively) and found 
that the evolution of the conductance distribution as a 
function of the length of the ID wire is governed by a 
Fokkcr-Planck equation. Thus the solution of such a diffu- 
sion equation gives the distribution of conductances p s {G) 
for a wire of length L, which is written as [11] 



Ps(G) 



G 2 



dy 



ye" 



x /cosh 2 / + 1 - 2/G 



where yo = arcosh(2/G — 1) and s = L/l. We make a 
couple of general comments on the scaling approach of 
Ref. [TU]. On the one hand, in order to obtain the above- 
mentioned diffusion equation, it is assumed that for a very 
small sample of length SL, the resistance average is fixed 
and depends linearly on SL. As a consequence the aver- 
age (—lnG) of the whole sample scales linearly with L: 
(— lnG) = L/l. On the other hand, one can see from eq. 
(UJ) that the distribution of the conductance is determined 
by the single parameter s. This is a remarkable fact since 
it means that all the information we need to give a statis- 
tical description of the conductance is the system length 
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Fig. 2: The average (— lnG) as a function of L for a wire with 
potential barriers whose separations and widths are generated 
accordingly to a Levy- type distribution with a — 1/2 (eq. © )• 
As we predict (— lnG) oc L 1 ^ 2 . The solid line is obtained from 
eq. (0. A good agreement is seen between numerics (squares) 
and theory. Inset: (—lnG) for a disorder distribution with 
a = 3/2. In agreement with our model, (—lnG) oc L, for 
a > 1, as in usual disorder systems. 



measured in units of the mean-free path, all other micro- 
scopic details arc irrelevant. This point was referred as a 
central-limit theorem for weakly disordered wires and eq. 
([1]) shows explicitly the single parameter dependence of 
the statistics of the transport problem. In relation to the 
previous comments, we shall see that for Levy-type disor- 
der (In G) does not scale linearly with L but still the con- 
ductance distribution is determined by the average (In G) 
and the exponent of the power-law decay of the disorder- 
configuration distribution. 

As a further motivation, let us illustrate the fact that 
the ensemble average (lnG) does not scale linearly with 
L in the presence of long-range disorder for a < 1. In 
fig. [5] we show the numerical results for (—lnG) for two 
different values of a. For the numerical simulations, we 
consider a succession of square-potentials barriers with 
fixed hight, but separations and widths randomly gener- 
ated from a Levy- type distribution jT^J- The statistics is 
collected across an ensemble of different disorder realiza- 
tions where the length of the system is fixed. For a = 1/2, 
we can see clearly from fig. [2] (main frame) that (— lnG) 
is not a linear function of L, whereas for a = 3/2 (inset of 
fig- 12) ( — lnG) oc L, as in the standard scaling approach. 
The solid line in fig. [5] is obtained from our theoretical 
model which we explain below. In order to describe the 
results shown in fig. [2] and calculate the conductance dis- 
tribution, it is essential to consider the strong fluctuations 
in the number of scatterers in a wire of length L due to 
the long tails of the density distribution of the disorder- 
configurations. For instance, if the tail of the density de- 
cays slowly as 

(2) 
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for large x and < a < 1, the first and second moment 
diverges. Actually, we shall concentrate precisely in these 
values of a as they exhibit the most interesting results for 
the conductance statistics. For the case 1 < a < 2, where 
the hrst moment is finite, the statistical properties of the 
conductance are similar to those described by the standard 
scaling theory (e.g., see the inset of fig. [5]). This situation 
can also be treated along similar lines of this work [13j . 

Let us start by introducing the probability density 
TIl(v) for the number of scattering- units v in a wire of 
length L. If the separations and widths of the scattering- 
units in the wire follow a law like eq. ([2} one can show, 
using the properties of the Levy distributions that 
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-q a , c (L/(2v) 1 ' a ), 
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for < a < 1, in the macroscopic limit L 3> c 1 /". Here 
q a)C is the probability density function (PDF) of the Levy 
distribution supported in the positive semiaxis that be- 
haves like p(x) for large values of x [see eq. ©]. The con- 
crete form of q a ^ c (x) is better expressed using the Fourier 
transform 

q a , c {k) = exp {-\k\ a (B9(k) + B*9(-k))) , (4) 

where 6 is the Heaviside step function and B = 
— cT(— o^e 4 " 2 ?, r being the Gamma function. Note that 
the result for IIl(i/), eq. depends only of the parame- 
ters c and a. All other details of the PDF p of the disorder 
configurations are washed out. This fact follows from the 
generalization of the central limit theorem to distributions 
with long tails (see footnote 1). 

In order to continue we introduce some notation: let 
(In G) v and (In G) l be the expectation values for samples 



1 Let us approximate the number of scatterers-units n by a contin- 
uos variable v such that the probability of having exactly n scatterers 

pn+ 1/2 

in a system of length L is given by J ^ il^(i/)dv. Therefore we 
have that 

pn+l/2 poo 

U L (u)du = / p 2n+ i(x)dx, 
10 J L 

where p m (x) is the mth auto-convolution of p(x) and represents the 
PDF of the total length after m drawings. 

Introducing the new variables 77 = L/(2u) 1 / a , x = x/(2n + l) 1 / a , 
and n L (77) = (a/2)L a V- a - 1 Ii{L a V-° ! /2). The previous relation 
reads 

poo poo 

U L (u)du = / (2n + l) 1 / a p 2n+1 ({2n + l) 1 ' a x)dx, 
A J A 

where A = L(2n + l)" 1 /". If we take now the limit of large L 
keeping A constant and use the generalization of the central limit 
theorem which states that 

lim m 1 ^ pm,{rn}/ a x) = q a ,c{x); 

m— >oo 

we obtain the relation 

poo poo 

lim / IIi(i7)dF= / q a c (x)dx, 
l ^°°Ja J A 

from which our result, Eq. (O, follows. 



with a fixed number of scatterers v and fixed length _L, 
respectively. From the standard scaling theory of local- 
ization (— lnG% is proportional to v. (— hxG) v = ais, a 
being a constant. Thus we have that 
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(- In G) L = / (-\nG),n L {u)du 
Jo 



(5) 
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where we have used eq. ([3]). Using the scaling property of 
the Levy distributions, namely c 1 / Q g Q , iC (c 1 / Q a:) = q a ^\(x) 
and making the change of variable z = L/(2ci>)~, we fi- 
nally obtain 

(-lnG) L = a — -/ z- a q aS {z)Az = a — I a , (7) 
c 2 J c 

where I a = (1/2) z~ a q a ^{z)Az. The last equality in 
i[7|) stands to emphasize the role played by the different 
elements (c, a) of the original distribution and the phys- 
ical characteristics (a, L) of the system. We have shown, 
therefore, that (— In G)l oc L a , in contrast to the linear 
behavior with L expected from the usual scaling theory. 
This is true for any value of < a < 1, while the standard 
linear behavior is recovered for a > 1. This was already 
illustrated in the inset of fig. [2] Moreover, In G for a < 1 
is not self-averaging. That is, in the limit of large samples 
we found that the ratio R = lrni£^oo var(ln G)/(ln G) 2 is 
non zero for a < 1, whereas R = for a > 1, as in the 
conventional scaling theory. 

We now calculate the conductance distribution. As we 
mentioned before, in the standard scaling theory the com- 
plete conductance statistics is fully determined (In G)l- 
We shall show that the same is true when the configura- 
tion of scattering-units follows a PDF with long tails, with 
the only additional information of the parameter a that 
characterizes the distribution tail. 

Assume that we have an ensemble of disordered wires 
with (— In G)l = £ whose scattering-units are distributed 
according to a long-tailed probability, as described before. 
The conductance distribution (G) is given by an integral 
that involves p s (G) 1 given in eq. ((T|) with the parameter s 
related to the number of scatterers, s = av, and in 
eq. ©, i.e., 



P S (G) 



PaAG)n L (u)du. 



(8) 



Now, using Eqs. J3]) and Q as well as the scaling prop- 
erties of the Levy distribution, we finally write the distri- 
bution of conductances as 



p d G ) = / Ps(a,t,z)(G)q a ,i(z)dz, 
Jo 



(9) 



where we have made again the change of variable z = 
L/(2cv)~ and defined s(a, £, z) = £,/(2z a I a ). We remark 
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that our result for P$(G) [eq. ©] depends only on two 
parameters (— lnG)z, = £ and a, as we have previously 
announced. Thus, other details of the disorder configura- 
tions are irrelevant. 

Next we show some specific examples of the distribution 
of conductances as given by eq. © ■ The theoretical re- 
sults are compared with numerical simulations using the 
model of potential barriers introduced previously. First 
we illustrate the fact that P$(G) is only determined by 
the values of £ and a. With this purpose we have calcu- 
lated numerically the conductance distributions from two 
differently generated ensembles of wires: the value of £ is 
the same in both cases, as well as the value of a, but the 
disorder-configurations of each ensemble follow different 
density distribution 0. We thus expect the same conduc- 
tance distribution for both cases. This is shown in fig. [3] 
where the distribution P^ (G) for each case (histograms in 
black and red) are statistically the same. In addition, in 
fig. [3J we plot in solid line our theoretical result, eq. ©• 
A good agreement with the numerical simulations is seen. 
On the other hand, we point out that the conductance 
distributions exhibit two peaks, at G ~ and G ~ 1. 
This is an unconventional behavior in ID disordered sys- 
tems and reveals the coexistence of insulating and metallic 
(ballistic) regimes. Interestingly, a similar behavior of the 
conductance distribution has been found in the so-called 
random- mass Dirac model |14j . where anomalously local- 
ized electronic states |15j might occur due to rare trapping 
disorder configurations. In the inset of fig. [3] we show a 
further example of the conductance distribution for wires 
with Levy- type disorder with a = 3/4. Again, the agree- 
ment between theory and numerics is very good. 

Finally, we show an example of P^{G) in the insulating 
regime, i.e., for very long wires. From the standard scaling 
theory, it is known that p s (G) follows a log-normal distri- 
bution whose variance is twice the mean value (— InG). 
An example of this standard case is shown in the inset of 
fig. 21 For long-range disorder-configurations the situa- 
tion is different. In the insulating regime, the integral in 
eq. © can be performed analytically using saddle point 
approximation. We plot the resulting distribution in fig. 
[4] (main frame) for a = 1/2 and the corresponding re- 
sult from the numerical simulation (histogram). We can 
see clearly that P^(lnG) is not Gaussian, as we have re- 
marked previously. Finally, we point out that in all previ- 
ous comparisons between theory and simulations there are 
no free-fitting parameters: given the value of the exponent 
a, all we need is (— lnG)z,(= £), which we extract from 
the numerical experiments. 

Conclusions.- Levy processes have been found in very 
different phenomena in nature, but recent transport exper- 
iments indicate the possibility of studying these processes 
in a controllable manner. Motivated by these experimen- 
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Fig. 3: (Color online) Conductance distributions from numer- 
ical simulations (histograms) and theory (solid line), eq. (0. 
Main frame: two histograms in black and red are shown (with 
different size bins for a better distinction). Each histogram is 
generated from different disorder distributions, pi(x) and p2{x) 
(see footnote 2), both with a = 1/2. In both cases £ = 1.1. 
P( (G) shows two peaks at G ~ and G ~ 1, revealing the co- 
existence of insulating and ballistic regimes. Inset: -Pj(G) for 
a — 3/4 with £ = 0.76 is shown. A good agreement between 
theory and numerics is seen in all cases. 



tal achievements, we have studied the statistical properties 
of the conductance across one-dimensional quantum wires 
whose disorder configurations are characterized by proba- 
bility densities with a power-law tail, i.e., Levy- type dis- 
order. We have shown that the conductance distribution 
is completely determined by only two parameters: the av- 
erage (In G) l and the scaling exponent a of the power-law 
tail. Therefore the distribution of conductances is univer- 
sal in the sense that all wires with the same (In G) l , but 
different disorder distributions obeying the same scaling 
exponent a have the same conductance statistics. Since 
we provide the complete conductance distribution other 
quantities of interest such as the moments of the conduc- 
tance and the shot-noise power can be calculated from 
our results. We have shown that for a < 1 the strong 
fluctuations of the conductance allow for the coexistence 
of localized and delocalized regimes and the fluctuations 
of In G do not exhibit self-averaging. For a > 1 the log- 
arithm of the conductance is a self-averaged quantity, as 
in the standard scaling theory. We think that the exper- 
imental observation of the different phenomena described 
here such as the universality of the conductance statistics, 
in optical materials or electronic transport experiments, 
should be of considerable interest. Finally, we leave for 
future investigations the extension of our results to higher- 
dimensional wires. 



2 Dcnsity distributions used for the numerical simulations: 
Pl(x) = 9i/2,cW = (c/v / 2^)x- 3 / 2 exp(-c 2 /2a;), p 2 (x) = 2~ 1 (1 + 
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Fig. 4: Conductance distribution for large disordered wires 
with a = 1/2 and £ = 226. We can see that P f (lnG) do 
not follow a Gaussian distribution, in contrast to the scaling 
theory prediction for the insulating regime, like the example in 
the inset. A good agreement is seen between analytical (solid 
line) and numerical (histogram) results. 



REFERENCES 

[1] Anderson P. W., Thouless D. J., Abrahams E. and 
Fisher D. S., Phys. Rev. B, 22 (1980) 3519. 

[2] Sims D. W., et. al., Nature, 45 (2008) 1098. 

[3] Mantegna R. N. and Stanley H. E., Nature, 376 
(1995) 46. 

[4] Levy P., Theorie de V addition des variables aleatoires, 
(Gauthiers-Villars, Paris) 1937; Gnedenko B. V. and 
KOLMOGOROV A. N., BookLimit distributions for sums 
of independent random variables (Addison- Wesley, Cam- 
bridge) 1954; Calvo I., Cuchi J. C, Esteve J. G. and 
Falceto F., J. Stat. Phys. , 141 (2010) 409. 

[5] Uciiaikin V. V. and V. M. Zolotarev V. M., Chance 
and Stability. Stable Distributions and their Applications 
(VSP, Utrecht, Netherlands, and references therein) 1999 

[6] Kohno H. and Yoshida H., Solid State Commun, 132 
(2004) 59 

[7] Barthelemy P., Bertolotti J. and Wiersma D. S., 

Nature, 453 (2008) 495. 
[8] Leadbeater M., Falko V. I. and Lambert C. J., Phys. 
Rev. Lett., 81 (1998) 1274; BOOSE D. and Luck J. M., 
J. Phys. A: Math. Theor., 40 (2007) 140405; Burioni 
R., Caniparoli L. and A. Vezzani, larXiv:1003.2161l 
Beenakker C. W. J., Groth C. W., and Akhmerov 
A. R., Phys. Rev. B, 79 (2009) 024204; A. Eisfeld, S. 
M. Vlaming S. M., Malyshev V. A. and Knoester 
J., Phys. Rev. Lett, 105 (2010) 137402 
[9] Mel'nikov V. I., Pis'ma Zh. Eksp. Teor. Fiz., 32 (1980) 
244 [JETP Lett, 32, 225, (1980)]; Dorokhov O. N.„ 
Pis'ma Zh. Eksp. Teor. Fiz., 36 (1982) 259 [JETP Lett, 
36, 318, (1980)]. 

[10] Mello P. A., J. Math. Phys., 27 (1986) 2876. 

[11] Beenakker C. W. J., Rev. Mod. Phys., 69 (1997) 731. 

[12] The incidence energy of electrons is five times the barriers 
height. For numerical generation of Levy random variables 
see Ref. [5] and Chambers J. M., Mallows C. L., and 



p-5 



